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Abstract. Slow dynamics of energy transfer between different phonon modes under the resonance conditions is considered. 
It may result in new effects in the inelastic and quasielastic neutron scattering spectra. 



INTRODUCTION 

Crystal lattice dynamics Ql 01 is a prototype theory for 
many-body physics as a whole. In this case all principal 
approximations and concepts can be justified accurately 
on the basis of adiabatic (Born-Oppenheimer) approxi- 
mation of quantum mechanics. Due to the smallness of 
typical atomic displacements in crystals u in compari- 
son with interatomic distances d one can pass rigorously 
from a problem of strongly interacting particles (atoms, 
ions, or molecules) to a problem of weakly interact- 
ing quasiparticles (phonons). In the leading order in the 
smallness parameter 77 = u/d the crystal lattice dynam- 
ics and thermodynamics can be described in terms of an 
ideal phonon gas (harmonic approximation) fll3]- With 
the temperature T increase the parameter 77 increases as 
well, however, due to a semiempirical Lindemann cri- 
terion (see, e.g., Q, y[ and Refs. therein) 77 ~ 0.1 at 
the melting point T = T m so higher-order {anharmonic) 
contributions to thermodynamic properties are usually 
small up to the melting temperature ||4J,q|]. This state- 
ment is true for the most of average characteristics. At 
the same time, for some peculiar modes anharmonic- 
ity can be crucially strong, especially in the vicinity of 
some structural transformations (see, e.g., computational 
results for Ba |01 and Zr 00], as well as inelastic neu- 
tron scattering data for A15 structure compounds |9] and 
for high-temperature bcc phases of Ti, Zr, and Hf lIlOlO . 
For these cases a standard phonon picture is not ade- 
quate, for instance, additional non-phonon peaks in the 
dynamical structural factor may appear |11]. A "slow" 
(low-frequency) lattice dynamics in the form of a "cen- 
tral peak" in quasielastic neutron scattering is also typical 
for these strongly anharmonic modes I9l fl2lfl3ll . 

Phonon picture can be broken for these modes since 
the corresponding effective potential V (u) turns out to 



be essentially nonparabolic 10, El E Q- In the 
vicinity of structural phase transitions (e.g., ferroelectric 
transitions 1 13] or martensitic transformations in metals 
la [Till ) this potential normally has several minima corre- 
sponding to several competing phases, the height of the 
barrier being much smaller than a typical cohesion en- 
ergy. It is not surprising that the harmonic approximation 
is an extremally poor starting point to describe this situ- 
ation. Here we discuss the following issue: is the phonon 
picture always adequate in a generic case when 77 << 1 
for all phonon modes? To answer this question it is in- 
structive to consider a crystal from the point of view of 
a contemporary theory of dynamic systems rather than 
from that of statistical thermodynamics 1 14]. Accor ding 
to the KAM (Kolmogorov-Arnold-Moser) theory fl5ll 
a system of noninteracting oscillators (that is, an ideal 
phonon gas) is stable with respect to small interphonon 
interactions (anharmonicity) except the cases of the res- 
onance where the ratio of phonon frequencies is close to 
the ratio of small enough integer numbers. One can ex- 
pect that this exceptional case is a situation where even 
small anharmonicities can lead to essentially nonpertur- 
bative effects in the lattice dynamics I14l[m[l7ll . 



RESONANCE PHENOMENA IN 
PHONON SUBSYSTEM 

The well-known Fermi resonance in the Raman molecu- 
lar spectra under the conditions of integer ratio of phonon 
frequencies I18l |2|] is a prototype example of these non- 
perturbative resonance phenomena. It was first discov- 
ered for CO2 molecule where the ratio of two molecular 
vibration frequencies oj\ and (Oi is close, by accident, to 
1:2. This leads to a quasidegeneracy of the states with 
energies h(£>\ and 2h(Oi. As a result the first harmonic of 



the first mode is strongly coupled with the second har- 
monic of the second mode which leads to the violation 
of symmetry-induced selection rules and appearance of 
additional lines in the Raman spectra. Fermi resonance 
is typical for many organic compounds, especially con- 
taining CH group, and is well investigated by molecular 
spectroscopists 1 19]. The theory of the Fermi resonance 
is rather simple since the ultraquantum limit T <C ho) 
(where 00 is a mean phonon frequency) is typical for all 
these cases and therefore there are only few states which 
are involved in the resonance. 

When we consider Fermi-resonance-like phenomena 
for crystals we immediately face several difficulties 
which make the problem much less trivial Jl4l fl6l fl7ll . 
First, the classical limit T > hod is interesting now and 
it turns out to be much more complicated than the ultra- 
quantum one since here we have an infinite number of 
relevant states. Second, now we have to satisfy the reso- 
nance conditions not only for the frequencies but also for 
the phonon wave vectors q. Then, there is a continuum of 
nonresonant phonons which is nevertheless also impor- 
tant working as a thermal bath and leading to the phonon 
damping and "Brownian motion" fluctuations. The most 
important point is the occurrence of essentially new phe- 
nomenon in the classical case, namely, a low-frequency 
dynamics due to energy transfer between the modes par- 
ticipating in the resonance; this effect is exponentially 
small at low temperatures (ultraquantum regime) when 
the corresponding modes are not excited. This "slow" 
dynamics corresponds to transitions between the states 
split by the anharmonic interaction under the Fermi res- 
onance conditions. In this sense it relates to the Fermi 
resonance like the electron paramagnetic resonance re- 
lates to the Zeeman effect. Experimentally it may result 
in the appearance of low-e nergy p eaks in the quasielactic 
neutron scattering spectra Ml^lla. lim . 

These new phenomena result from a correlated char- 
acter of atomic displacements corresponding to different 
phonon branches under the Fermi resonance conditions. 
In a generic case the smallness of adiabatic parameter 
justifies the applicability of the perturbation theory for 
the consideration of anharmonic effects l2flH 3l. From the 
technical point of view, it means the decoupling of the 
higher-order phonon correlators according to the Wick's 
theorem, or, in classical terms, we suppose that the ini- 
tial phases of the interacting phonons are random. This 
is a standard approximation which is used to calculate 
the temperature d epen dence of phonon frequencies and 
dampings 12ll 122112311 . However, when the phonon fre- 
quency ratio turns out to be integer (and provided that 
some symmetry conditions are satisfied) a well-known 
nonlinear effect, the "phase locking", appears l24ll and 
the perturbation theory fails. In particular, it leads to an 
essential enhancement of multiphonon contributions to 
the dynamical structure factor in comparison with the 



perturbative treatment, and to the appearance of qua- 
sistatic atomic displacements. This is a consequence of 
subtle phase coherence phenomena which are completely 
igonored by a standard perturbative treatment. The situa- 
tion is similar to the Anderson localization in disordered 
systems which also cannot be described by any consider- 
ation of the average Green function |25]. 

To illustrate general mechanisms consider first the 
simplest case of the resonance, namely, phonons with 
the wave vector qo = (3,5,5) and the frequency ratio 
1:2 (HI 0- It corresponds to BCC phases of alkali and 
alkaline earth metals (for instance, for potassium the fre- 
quencies of longitudinal and transverse phonons in this 
point at T = 4.2K are 0.27 ±0.01 and 0.55 ±0.01 of ionic 
plasma frequency, respectively). Transverse phonons are 
double degenerate, however, microscopic estimations of 
the corresponding anharmonic coupling constants fl^l 
show that only one transverse branch participates in the 
resonance since the coupling constant for another one is 
four orders of magnitude smaller. So we can consider 
a simplified model with the interaction of one longitu- 
dinal branch (with the displacement field u(r,t)) and 
one transverse one (with the displacement field v(r,f)). 
Taking into accoun only resonant interaction term V = 
X J druv 2 one can write a set of equations of motion, 

f u + (0 2 {-N)u + 2y{-iV)u + lv 2 = 

\ v + i2 2 (-/V)v + 2r(-/V)v + 2AMv = '" ! 

where 7, Y are corresponding phonon damping parame- 
ters. We try the solutions of Eq.Q in the form 

' m — A(r,t) exp/ [qor — 00 (qo) t] + 
B(r,f)exp/[q r+o)(qo)f]+c.c. 
v = C (r, t ) exp/ [q r + ffl (q ) t /2] + 

k D (r,t) exp/ [q r - o(q )f/2] +c.c. 

Amplitudes A,B,C,D are slowly varying in space and 
time functions. It is important that the wave vector qo 
is equivalent to — 2qo since exp(3/qor) =1 for the crys- 
tal lattice sites. Substituting Eq.(|2} into Eq.Q, taking 
into account only resonant terms and only leading ap- 
proximations in the anharmonic smallness parameter (for 
more details, see Ref. 1 1711 ') 

dt \d(l ) q 4fflo \dq a dqp ) dx a dxp 

+y A + /AC* 2 = 
* dc_i(d®.) vc-i L ( d2 °> 2 ) d2c ^ 

dt 2 \ dq J 8fflo \dqadqp ) q dxadxp 

+ (r + i -f)c- 4/AA*C* = 

where co = 2D. + v (V << 03) ,A = A/2«o and subscript 
"0" means q = qo. Equations for B and D differ from 
Eq.Q by the replacement v — > — V, A — > —A. 

One can demonstrate 1 1 711 that at small enough v (re- 
alistic estimations for alkali metals show that "small" 



means 5-7% of CO) purely sine waves A,B,C,D = const 
turn out to be unstable with respect to a self-modulation 
and the solitons of the envelopes can form. This mod- 
ulation is connected with a slow (in comparison with a 
characteristic phonon times) dynamics of energy trans- 
fer between longitudinal and transverse phonons. Nu- 
merical simulations of the thermal noise effects 11261 12711 
show that the latter do not suppress this slow dynamics. 
It appeared that there are two limit circles in this system 
(which correspond to two phase locking regimes with 
different relative phases for A and C waves) and this en- 
ergy transfer dynamics can be described as a stochastic 
resonance between these two limit circles i27ll . 

Similar phenomena can be also considered for acous- 
tic (long-wavelength) phonons as an energy transfer be- 
tween two ultrasound waves with different polarization 
vectors and integer ratio of sound velocities 12811 . Al- 
kali metals near the melting point (1:3 ratio for differ- 
ent transverse sound waves propagating into <110> di- 
rection) or Wi_;tRe x , Moi_ x Re v alloys (1:1 ratio for the 
same sound waves) might be an interesting examples. 
Computer simulations 12811 demonstrate that, depending 
on the initial phases of the waves, this energy transfer 
may be both chaotic and quasiperiodic. It would be inter- 
esting to check this prediction experimentally. 

It is worthwhile to note that for the case of 1:2 fre- 
quency ratio considered above the resonance conditions 
for phonon wave vectors can be satisfied only in some pe- 
culiar points of the Brillouin zone. For the case of 1 :3 res- 
onance they can be satisfied in a generic case; for high- 
symmetry directions the number of waves participating 
in the resonance coupling can be very large. FCC La is 
an interesting example of such resonance (below we fol- 
low our work |29]). The lattice dynamics of FCC La is 
characterized by a drastic nonmonotonicity of the dis- 
persion curves ffl(q) for transverse phonons in q || (111) 
direction and by a significant temperature dependence 
of their frequencies, i.e. by strong anharmonic interac- 
tions 1 30]. As it was discussed in Ref. 113 Oil this behavior 
can be accounted for by the electron-phonon interaction 
(Kohn anomalies). One can see from these experimental 
data that the ratio of longitudinal to transverse phonon 
frequencies in (111) direction varies around the integer 

3, the value of wave vector qo varying 



£»i(qo) 



ratio 

with the temperature within sufficiently wide limits. This 
makes it possible to shift the resonance watching simul- 
taneously the change in the lattice dynamics behavior. 

Let us set up the simplest model to describe the res- 
onance effects under consideration in the lattice dynam- 
ics of FCC La. We shall proceed with the equations of 
motion for the amplitudes of longitudinal (u) and trans- 
verse (v) phonons allowing only for the "resonance" an- 
harmonic interaction V = X f druv^: 



oh 

• qV 



u + co u 



-Av 3 = 
-3Xuv 2 = 



(4) 



Because of the FCC lattice symmetry we should 
consider phonons propagating in four equivalent direc- 
tions (111), (111), (III), (111); the relevant wave vec- 
tors meeting condition ffl(q) = 3£2(q) will be denoted as 
q ; (j = 0, 1,2,3). Similarly to Ref. 1 17] we try the solu- 
tions of the equations (0) as 

3 

u = Y. (^i ex P'(qi r _ cot)+Bje\pi(qjr+ cot)) +c.c, 

7=0 
3 

v= £ (C/expz(q y r — Qf) + D/expi(q/r + Qf)) +c.c. 

(5) 

where Aj,Bj,Cj,Dj are slowly varying (due to X small- 
ness) functions of r and t, i.e. the envelopes of the 
phonons under consideration. Substituting Eqs. Q into 
Eqs. © and neglecting the second time derivatives of 
the envelopes as well as the nonresonance terms we ob- 
tain the following set of equations 



A + iA(6D*D* 2 D* 

B - zA(6CfC|C| - 
C + 9/A(2fi|qC 3 * 



3C 2 D*- 







6C I CjD* 

.7=1 

3£>2q+6A, I Djq) = o 

7=1 



-2B* 2 C\C\ 



■2B\C\C\- 



2Aq L CjDj+2Do L A jCf 



DiB : 



0"0" 



Do - 9iA(2A*D* 2 D 3 



2D ZDjB*j) 

7=1 



2A* 2 D\Dl 



2A\D\D* 2 - 



2B I D* j Cj + 2C Z BjD* r 



C 2 A* 



■2C I CjA*) 

7=1 



0; 



(6) 

the rest twelve equations are obtained from (|6jl by cyclic 
permutations of the indices. Here A = 

The dynamics of the energy transfer in this system has 
been investigated in Ref. 1 29] by numerical simulations. 
It was shown that depending on the initial phases of the 
involved phonons this energy transfer can be either reg- 
ular or chaotic. One can expect from the dimension con- 
siderations that a typical frequency of this slow dynamics 
will be of order of 



co 



(7) 



Here M is the ion mass, a is the lattice constant, u 2 is the 
average square of atomic displacements. For the room 
temperature CO* w 10~ 3 O> can be taken for the estimation. 
However, it was shown that in reality it is much larger, 



A(0 ~ 10 2 (O* , due to a large number of interacting waves 
(32 real fields). 

Similarly to Ref. 11711 (see also above) it can be shown 
that the account for intermode or non-resonant intramode 
anharmonicities may result in the appearance of the 
peaks in quasielastic neutron scattering spectra with a 
frequency width about Ao. The largest contributions re- 
sult from the lowest-order (three-phonon) anharmonic 
processes. For example, the term V = fif afrv 3 in the 
potential energy leads to the appearance of the "low- 
frequency" contribution to the transverse phonon field, 
Sv = —3jJ.v 2 /Q. 2 . Assuming all the amplitudes Aj and Cj 
in the equation (0 to be constant we would obtain static 
contributions to v with the wave vectors Q a = q, ± qj + g 
where g is a reciprocal lattice vector (which can be, in 
particular, equal to zero). One can expect therefore the 
peaks in the quasielastic neutron scattering spectra near 
the wavevectors Q ( , which reflect the dynamics of the 
envelopes. In a special case considered above for the fre- 
quency ratio 2:1 one of these vectors coincides with the 
vector q,- but for the ratio 3:1 it is impossible. The vectors 
Q„ for i = j are just the reciprocal lattice vectors but the 
rest of them should be temperature dependent due to the 
temperature dependence of the phonon frequencies shift- 
ing the Fermi resonance point. It can be important for the 
experimental verification of the effects under considera- 
tion. A similar contribution of order of u 2 appears in the 
longitudinal phonon field because of the potential energy 
term V = V f dru 3 . 

CONCLUSIONS 

To conclude, we predict a new class of nonperturbative 
anharmonic phenomena in lattice dynamics under reso- 
nance conditions in the phonon system. This resonance 
results in the instability of phonons and complicated pic- 
ture of energy transfer between the modes in the reso- 
nance. It can be experimentally investigated by inelastic 
and quasielastic neutron scattering as well as by acoustic 
methods. There is an interesting open issue what is the 
role of these phenomena in the development of lattice 
instabilities of the crystals, in particular, in the melting 
0. This question is motivated by two general remarks: 
(i) anharmonic effects in a generic case are small up to 
the melting point and (ii) instability of a generic weakly 
anharmonic dynamical system, according to the KAM 
theory, is connected with the resonances. Therefore it 
might appear that the phenomenon under consideration 
is not a kind of exotic but important for any crystal. 
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